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Chapter  1 

The  Random  Structural  Eigenvalue  Problem 


1.1  Introduction 

In  this  chapter  we  will  take  the  first  step  in  the  direction  of  incorporating  in  the  d}rnamic  analysis 
of  the  structures  the  uncertainties  associated  with  the  values  of  the  connection  stiffnesses. 

Under  the  assumption  that  the  moment-rotation  relationship  for  the  joints  is  linear,  each 
connection  stiffness  is  represented  by  a  single  parameter,  i.e.,  the  slope  of  the  M  —  0  curve.  This 
parameter  should  be  determined  from  a  series  of  tests  performed  for  each  type  of  connection  and 
for  a  given  size  of  the  structural  elements  that  converge  to  the  joint.  The  results  presented  in 
[1],  [2]  and  [3]  may  be  mentioned  as  examples  of  such  experimental  investigations.  In  particular, 
Altman  et  al.  [3]  identified  the  most  important  parameters  that  control  the  initial  stiffness  in 
bolted  beam-column  connections.  Any  change  in  the  value  of  these  parameters  will  cause  a 
variation  in  the  linear  relaticmship  assumed. 

Although  using  the  test  results  it  is  possible  to  assign  a  value  to  the  rotational  stiffness  of 
a  given  connection  configuration,  one  cannot  guarantee  that  this  value  will  be  the  same  for  all 
the  connections  of  the  same  type  in  a  given  structure.  The  values  of  the  connection  stiffness 
obtained  from  a  test  shoul^be  regarded  as  the  mean  value  of  the  stiffness  of  connections  with 
the  same  configuration.  There  will  always  be  a  degree  of  dispersion  due  to  the  variations  in  the 
several  parameters  that  control  the  connection  behavior.  Therefore,  it  is  logical  to  consider  the 
connection  stiffness  as  a  random  variable  in  the  structural  model.  In  this  case  both  the  stiffness 
and  the  mass  matrix  will  have  random  coefficients,  since  both  are  function  of  the  fixity  factors, 
which  in  turn  are  defined  in  terms  of  the  connection  stiffness.  The  finite  element  equations  of 
motion  of  the  structure  become  a  set  of  ordinary  differential  equations  with  time-independent 
random  coefficients.  Although  the  equations  of  motion  are  random  they  are  still  linear  and  hence 
they  could,  at  least  in  principle,  be  decoupled  and  solved  by  modal  analysis.  The  eigenvalue 
problem  associated  with  the  equations  of  motion  also  becomes  random,  and  the  eigenvalues  and 
eigenvectors  of  the  system  end  up  being  random  variables  as  well.  The  object  of  this  chapter 
is  to  study  this  random  algebraic  eigenvalue  problem  associated  with  a  structural  system  with 
random  flexible  joints.  The  dynamic  response  of  these  systems  will  be  addressed  in  the  following 
chapter. 
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1.2  Connection  Stiffness  as  a  Random  Variable 


The  determination  of  the  probability  associated  with  a  particular  value  of  stiffness  requires  the 
knowledge  of  either  its  distribution  function  or  its  density  function.  In  general,  the  choice  of 
a  distribution  to  describe  the  probabilistic  behavior  of  a  physical  variable  is  governed  by  the 
nature  of  the  phenomenon.  In  the  particular  case  of  a  variable  that  represents  the  total  effect  of 
several  random  causes,  the  central  limit  theorem  leads  to  the  conclusion  that  its  distribution  is 
asymptotically  normally  distributed.  As  it  was  mentioned  before,  there  are  a  number  of  factors 
that  contribute  to  the  variability  of  the  stiffness  of  a  given  connection.  Consequently,  it  is 
reasonable  to  assume  that  the  stiffness  values  can  be  modelled  as  random  variables  with  normal 
or  Gaussian  distribution. 

In  our  case,  the  value  of  the  connection  stiffness  can  be  regarded  as  the  sum  of  a  deterministic 
component  and  a  component  representing  random  perturbations  due  to  the  variability  in  the 
controlling  parameters.  These  variations  can  be  caused,  for  example,  by  fluctuations  in  the 
material  characteristics,  the  geometry  of  the  connection,  the  assembling  method  or  the  quality 
control  of  the  process.  Therefore,  the  stiffness  value  of  the  i—th  structural  joint  can  be  represented 
as  follows: 


Ki  =  Ki  +  oi  (1.1) 

where  Ki  is  the  random  rotatimial  stiffness  of  the  joint,  ICi  is  a  deterministic  constant  and  Qi  is  a 
random  perturbation.  In  order  to  simplify  the  analysis  the  mean  value  of  the  random  perturbation 
will  be  incorporated  into  the  deterministic  component  so  that  a,-  can  be  regarded  as  a  zero-mean 
random  variable.  Moreover,  if  we  assume  that  the  probability  distributitm  of  the  connection 
stiffiiess  is  normal,  then  we  conclude  that  Oi  is  also  normal,  since  normal  random  variables 
remain  normal  under  linear  transformations.  With  the  above  considerations  it  is  straightforward 
to  show  that  the  first  two  moments  of  the  random  variables  a,-  and  Ki  are  given  by: 

^  =  0  ;  E{Ki}  =  Ki  (1.2) 

£{(.()’}  =  <  =  <^,  ;  =  ■^,  +  a?  (1-3) 

Hence,  the  mean  value  of  the  stiflhess  is  Ki ,  and  its  variance  a^Ki  is  equal  to  the  variance  of 
the  random  perturbation  .  These  two  quantities  define  completely  the  probability  distribution. 


1.3  Second  Order  Perturbation  Technique 

In  the  formulation  of  a  probabilistic  finite  element  method  based  on  the  second  order  perturbation 
technique,  each  random  variable  is  expanded  about  its  mean  value  amd  terms  up  to  second 
order  are  retained.  The  rates  of  change  of  the  eigenproperties  with  respect  to  the  fixity  factors 
calculated  in  the  previous  chapter  will  be  used  to  obtain  expressions  for  the  mean  and  variances 
of  the  eigenproperties  in  terms  of  the  first  and  second  moments  of  the  random  stiffnesses.  For 
this  reason,  this  type  of  approach  is  called  a  second  —  moment  analysis.  The  inherent  limitation 
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of  this  formulatiim  is  that  the  statistical  variations  of  the  random  variables  have  to  be  small  in 
order  to  obtain  acceptable  accuracy,  as  it  is  the  case  with  all  perturbation-based  methods. 

If  the  random  parameters  of  the  structure  are  substituted  by  their  means  or  expected  values, 
we  obtain  an  averaged  version  of  the  structural  system.  Therefore,  using  the  mean  values  of  the 
joint  stiffness  to  define  the  system  matrices,  the  solution  of  the  associated  eigenvalue  problem 
results  in  a  deterministic  set  of  eigenproperties.  It  will  be  shown  later  that  these  eigenproperties 
coincide  with  the  mean  values  of  those  obtained  through  a  perturbation  method  based  on  a  first 
order  expansion.  Hence,  to  study  the  difference  between  the  deterministic  eigenproperties  of 
the  averaged  system  and  the  mean  values  of  the  eigenvalues  and  eigenvectors  associated  to  the 
random  system,  it  is  necessary  to  include  at  least  second  order  terms  in  the  expansions. 

The  eigenvalue  problem  for  a  structural  system  with  N  unconstrmned  degrees  of  freedom  and 
R  random  flexible  joints  is: 


{[KiK)]  -  XiiK)  [M(if)]}  =  0  ;  i  =  1,2, ...,JV  (1.4) 

where  [if]  and  [Af]  are  the  stiffness  and  mass  matrices  whose  coefficients  are  function  of  the 
random  variables  ifi,  A'a,...,#^^.  These  variables  are  expressed  in  vector  form  as  follows: 

K  =  {KxK2...KRf  (1.5) 

The  eigenvalues  Aj  and  eigenvectors  ^  become  nonlinear  functions  of  the  variables  K\,  —iKr 

and  hence  they  are  also  random  quantities.  The  spedflc  form  of  the  nonlinear  functions  cannot  be 
determined,  except  for  trivial  cases.  Nevertheless,  assuming  that  the  variables  Ki  are  constrained 
to  small  fluctuations  about  their  mean  values,  we  can  express  the  eigenproperties  as  Taylor  series 
expansions  in  terms  of  the  random  stifbess  parameters.  The  coefficients  of  the  expansions  are 
evaluated  at  the  mean  value  vector  defined  as  follows: 


k  =  {kxK2...KRf 

Therefore,  the  eigenvalues  cSn  be  approximated  as: 


\i(K)  =  A,(i)  +  E  ^ 


R  R 


+  5EE 

msl  nsl 


dKr^dKn 


K=k 


K=k 

{Km  -  km){Kn  -  Kn) 


Introducing  the  notation: 


A,  =  A,(ir) 

dKmlK^k 


(1.6) 


(1.7) 


(1.8) 

(1.9) 
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(1.10) 


A."  = 

•mn 


dKrr,  dKn 


and  expressing  the  random  perturbations  from  the  mean  values  of  the  stiffness  coefficients  as: 


—  Ktn  (1.11) 

equation  (  1.7)  can  be  written  as: 

R  R  R 

Af  =  Aj  +  ^  A/^  Om  +  2  ^  ^  (1-12) 

m=:l  msl  n=l 

Proceeding  in  a  similar  way,  the  second  order  Taylor  series  expansion  for  the  eigenvectors  can  be 
written  as: 


4i  =  ^.  +  E  +  E  E  “m  “n 

msl 


R  R 

CE 

msl  nsl 


(1.13) 


where: 


-  ^i(k) 


(1.14) 


=  ^. 


- 


OKmdKn 


K^k 


(1.15) 

(1.16) 


In  the  previous  chapter  we  defined  the  Taylor  series  expansion  for  the  eigenvalues  and  eigenvectors 
in  terms  of  the  so-called  fixity  factors  /<,■  .  These  non-dimensional  factors  allowed  us  to  obtain 
compact  expressions  for  thd^coefficients  involved  in  the  calculation  of  the  rates  of  change  of  the 
eigenproperties.  At  first  sight,  it  seems  convenient  here  to  make  use  of  the  same  expansions. 
However,  recalling  that  the  fixity  factors  are  defined  in  terms  of  the  joint  stiffness  as: 


-  =  TTm, 

we  observe  that  if  the  variables  Ki  are  considered  to  be  normal  random  variables,  the  same 
assumption  cannot  be  extended  to  the  factors  m  ,  since  equation  (  1.17)  does  not  define  a 
Himr  tnaffonnatioii.  Therefore,  if  we  are  interested  in  maintaining  the  Gaussian  distribution 
assumption  to  take  advantage  of  its  properties,  the  expansions  have  to  be  written  in  terms  of  the 
variables  .  However,  we  can  still  employ  the  available  expressions  for  the  rates  of  change  of 
the  eigenproperties  making  use  of  the  chain  rule: 


;  /  _ 


\I 

dKm 


(1.18) 
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(1.19) 


_  xl 

dK^  ~  aiir„ 


A/i 


d^A. 


dKradKn 


=  VL  (te)"  +  S®  for  m  =  n 
=  1^  si?"  for  m  n 

\  •mm  vivfvt  v/vfi  ' 


(1.20) 


,// 


^  = 

^«mfi 


d^Ai 


dKmdKn 


in  which: 


iff” = 

(Jir„  +  3^i)’ 

a»<.„  _  -«¥ 

+ 


(1.21) 

(1.22) 

(1.23) 


To  obtain  the  expected  valnes  of  the  i-th  eigenproperties,  it  is  necessary  to  apply  the  expected 
value  operator  to  the  expansions  (  1.8- 1.11).  Considering  equation  (  1.2),  we  can  write: 

R  R 

E{Ai}  =  A.  +  5  z  I;a£L  £<<>..0.)  (1.24) 

W)  =  *  +  5  E  Eiii',  £{a«o.}  (1.25) 

^  msl  nsl 

FVom  these  equations,  we  can  see  that  the  expected  values  differ  from  the  zero-order  terms 
that  are  solution  of  the  averaged  problem,  and  the  difference  is  a  linear  function  of  the  covariances 
of  the  random  variables. 

If  the  random  variables  are  considered  to  be  uncorrelated,  then  it  is  possible  to  write: 


E{an,a„}  =  Srnn<Ti  (1.26) 

where  Smn  is  the  Kroenecker  delta,  and  <t\  is  the  variance  of  the  random  variable  >  With 
this  assumption,  equations  (  1.22-  1.23)  take  the  form: 

e{Ai}=S.  +  i  EA?;.*'  (1.27) 

^  m«l 

e{*i)  =  ^.  +  \'t*L<  (1.2*) 
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The  variance  of  the  t  -  th  eigenvalue  is  given  by: 

=  f;{(A*  -  E{Xi}f}  (1.29) 

Substituting  the  corresponding  expressions  for  Xi  and  given  by  equations  (  1.12)  and 

(  1.4),  respectively,  and  considering  the  linearity  of  the  expected  value  operator,  we  can  write: 

Imsl  rsl  asl  ) 

vms:l  f=l  *=1  } 

{R  R  R  H  ^ 

E  (1-30) 

matl  natl  r=l  **1  ) 

If  we  consider  the  properties  of  normal  zero-mean  random  variables,  we  are  lead  to  the  conclusion 
that  the  fourth  and  fifth  terms  in  the  above  expression  must  vanish,  since  they  involve  odd 
moments: 

*  £  { Otnt  Otr  }  —  0  (1*31) 


E{a„,E{ara,})  =  £{0^}  E{arQ,}  =  0  (1.32) 

and  equation  (  1.28)  reduces  to: 

m=l  nsl 

E  E  E E + 

msl  nasi  rsl  «sl 

R  R  R  R 

+  -  E  E  E  E  ^  }  (1.33) 

mal  nal  rsl 
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The  joint  moment  of  four  normal  random  variables  can  be  written  in  terms  of  lower  order  moments 
as  follows: 


E  {amOtnQrQ*}  =  {  <»«  }  f?  {  O,  Otj  }  +  £  {  Om  }  £  {  On  O,  }  + 


+  £  {  Om  Qt«}  £  {  On  Otr  }  (1-34) 

and  the  multiple  summati<Mi  in  the  last  term  of  equation  (  1.30)  can  be  expressed  in  the  following 
form: 

R  R  R  R 

£  H  E{  am  an  Or  a,}  = 

m=l  nsl  rsl  «sl 

msl  nsl  rsl  csi 

+  \"^)  (1-35) 


Substituting  equation  (  1.33)  into  equation  (  1.30),  we  finally  obtain: 

R  R 

<  =  E  Z 

,  R  R  R  R 

+  iEEEE 

msl  nsl  rsl  asl 

The  expression  for  the  eigeiyector  variance  can  be  obtained  following  an  analogous  derivation: 

‘"Ji  =  E  E  ^;{0tman}  + 

msl  nsl 

+  j  E  E  E  E  +  ^".C)  £«“-«.)  (1-37) 

msl  nsl  rsl  jsl  '  ' 

If  we  assume  that  the  random  variables  are  uncorrelated,  then  equations  (  1.34)  and  (  1.35) 
reduce  to: 


(^imr  4"  ^imt  ^iir  )  {  ®»n  On  }  ^  {  Or  Ol#  }  (1.36) 


=  E  (^)’  j  E  E  <  <  (1-38) 

msl  msl  nsl 

=  E  (*1)’  oil.  +  5  E  E  (♦!)’  ’i  (1-59) 

msl  '  '  ^  msl  nsl  '  ' 
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1.4  Numerical  Examples 

In  the  first  part  of  this  section  we  will  validate  the  expressions  for  the  first  and  second  order 
moments  of  the  eigenproperties  derived  in  this  chapter.  For  this  we  will  consider  the  plane  frame 
illustrated  in  Figure  (2.4)  and  we  will  compare  the  values  of  the  moments  obtained  with  the 
2nd  order  perturbation  expansion  method  against  the  results  of  a  Monte-Carlo  simulation.  The 
Monte-  Carlo  simulation  (hereafter  referred  to  as  MCS)  is  a  very  powerful  method  to  determine 
the  performance  of  systems  with  random  parameters.  The  method  consists  in  the  generation  of 
a  set  of  systems  derived  from  the  original  system  by  assigning  values  to  the  random  parameters. 
All  the  systems  are  later  statistically  processed  using  techniques  of  sampling  and  parameter 
estimation  [4,5].  The  values  assigned  to  the  parameters  of  the  systems  are  obtained  by  means  of 
an  algmithm  that  generates  random  numbers.  This  requires  to  select  the  probability  distribution 
function  that  governs  the  behavior  of  these  parameters. 

Once  the  accuracy  of  the  expressions  obtained  with  the  Second  Order  Perturbation  Method 
(from  now  on  designated  as  SOPM)  is  established,  they  will  be  used  to  obtain  the  probability 
density  function  of  the  eigenvalues  and  elements  of  the  eigenvector.  This  will  allow  ns  to  assess  the 
level  of  dispersion  in  the  dynamic  properties  of  the  structure  introduced  by  the  random  variation 
of  the  sti&ess  coefficients.  The  two-level  frame  in  Figure  (2.16)  will  be  used  for  this  purpose. 

1.4.1  Example  No.  1:  One-Story  FVame 

The  subroutine  GGNSM  of  the  STAT-IMSL  libtaiy  [6]  was  used  to  generate  two  sets  of  random 
numbers.  Both  sets  of  numbers  have  a  covariance  matrix  equal  to  the  identity  matrix.  A  pair 
of  50  uncorrelated  standard  normal  random  deviates  -  1,2,. ..,50)  with  zero  mean 

value  and  variance  equal  to  one  were  generated  in  this  way.  It  is  well  known  that  normal  random 
variables  remain  normal  under  linear  transformations.  Therefore,  these  variables  can  be  used 
to  generate  another  set  of  random  normal  variables  =  1,2,..., 50)  with  mean  value 

equal  to  k  and  standard  deviation  tri, .  The  rdationship  between  the  two  sets  is: 

+  fc  =  (Pfcili  +  l)i  ;  »  =  1,2,...,50  (1.40) 

k2i  =  <rkQi  +  k  =  {pkQi  +  1)*  ;  t  =  1,2,...,50  (1.41) 

Fifty  pairs  of  numbers  representing  the  stiffness  coefficients  of  the  rotational  springs  at  the 
beam’s  ends  were  generated  in  this  way.  The  mean  value  of  the  coefficients  k\^  ,  k^^  is  such  that 
the  corresponding  fixity  factor  is  0.50.  Two  different  sets  of  pair  of  coefficients  were  generated 
for  two  values  of  the  coefficient  of  variation  (c.o.v.)  of  the  stiffiiess  parameters,  pk  =  0.10  and  pk 
=  0.20. 

Tables  (1.1-6)  present  the  mean  values  and  standard  deviations  for  the  eigenvalues  and  eigen¬ 
vectors  of  the  first  two  modes  of  the  frame  in  Figure  (1.1)  calculated  with  the  two  approaches. 
The  frame  was  modelled  with  only  three  elements  and  6  dof  because  the  aim  of  the  analysis  is  to 
compare  the  accuracy  of  the  perturbation  method  without  bring  concerned  about  the  accuracy 
of  the  calculated  rigenproperties.  Tables  (1.1)  and  (1.2)  show  the  expected  value  and  standard 
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deviaticMi  of  the  lower  eigenvalue  of  the  frame,  respectively.  The  results  shown  under  the  Monte- 
Carlo  column  are  the  sample  mean  and  sample  standard  deviation.  All  the  tables  include  the 
limiting  values  corresponding  to  a  confidence  interval  of  95%.  It  can  be  seen  that  the  results 
obtmned  via  the  SOPM  always  lie  within  the  confidence  interval.  Note  also  that  according  to  the 
table,  doubling  the  standard  deviations  of  the  connection  stiffiiess  produces  the  same  effect  on 
the  standard  deviations  of  the  eigenvalues.  The  same  phenomenum  can  be  observed  in  the  c.o.v. 
of  the  eigenvalues  in  Table  (1.3).  This  reveals  an  almost  linear  probabilistic  dependence  between 
the  stiffness  of  the  connections  and  the  eigenvalues  of  the  structural  system.  Ebcamining  Table 
(1.3)  it  becomes  apparent  that  the  structural  system  filters  out  the  effect  that  the  uncertainties 
in  the  stiffness  of  the  connections  have  on  the  dgenvalnes  of  the  system.  For  instance,  when  the 
c.o.v.  of  the  stiffness  of  the  joints  is  0.10,  the  c.o.v.  of  the  first  and  second  eigenvalue  is  0.0139 
and  0.0378,  respectively. 

Tables  (1.4-6)  show  the  comparison  in  the  statistics  of  the  first  eigenvector  calculated  using 
the  SOPM  and  the  MCS  technique  for  two  values  of  the  c.o.v.  of  the  connection  stiffiiess.  It  is 
illustrative  to  discuss  some  characteristics  of  the  results  presented  in  these  tables.  The  expression 
for  the  expected  value  of  the  t  -  th  eigenvector  ^ven  by  equation  (  1.28),  includes  the  2nd  order 
derivatives  of  the  eigenvector.  It  is  shown  in  Ch^ter  3  of  ref.[7]  that  these  2nd  order  derivatives 
are  obtained  in  terms  of  a  linear  combination  of  the  eigenvectors  of  the  deterministic  system 
associated  with  the  mean  values  of  the  stiffiiess  co^cients.  Due  to  the  symmetric  configuration 
of  the  structure  under  consideration,  these  deterministic  eigenvectors  (or  modal  shapes)  are  dther 
symmetric  or  antisymmetric.  Thus,  the  expected  values  of  the  dgenvectors  of  the  random  system 
are  also  either  symmetric  or  antisymmetric.  The  same  sitnation  occurs  with  the  variance  of  the 
eigenvectors.  The  expression  for  the  variance  of  the  t  -  th  agenvector  is  given  by  equation 
(  1.39),  and  it  includes  the  squares  of  the  1st  order  derivatives  and  the  squares  of  the  2nd 
order  derivatives.  Therefore,  for  this  particular  mcample  the  variance  vectors  are  also  symmetric. 
On  the  other  hand,  the  results  of  MCS  are  obtained  from  the  statistical  processing  of  a  set  of 
structures  which  are  not  necessarily  symmetric,  since  the  stiffiiess  coefficients  are  uncorrelated 
random  variables.  Therefore,  for  a  sample  of  finite  size  the  expected  values  and  variances  of  the 
eigenvectors  do  not  necessagily  show  a  strict  symmetry.  In  spite  of  this  fact,  the  mean  values 
and  standard  deviations  of  all  the  elements  of  the  eigenvector  calculated  using  the  SOPM  fall 
inside  the  95%  confidence  interval.  The  c.o.v.  for  the  first  mode  of  vibration  are  presented  in 
Table  (4.6).  The  table  shows  that  the  modal  dof  associated  with  the  horizontal  displacement  of 
the  ends  of  the  beam  (1st  and  4th  elements)  ue  practically  insensitive  to  the  uncertainties  in  the 
stiffiiess  of  the  connections.  Moreover,  there  is  an  almost  linear  statistical  relation  between  the 
components  of  the  rigenvector  and  the  stiffiiess  of  the  joints.  In  this  as  well  as  in  the  previous 
tables  the  agreement  between  the  perturbation-based  results  and  the  Monte-Carlo  technique  is 
remarkable. 

1.4.2  Example  No.  2:  Two-Story  Plane  FVame 

The  next  example  is  the  two-story  plane  frame  modelled  with  42  dof  shown  in  Figure  1.2.  The 
stiffiiess  coefficients  of  the  connections  at  the  ends  of  the  beams  at  the  first  level,  denoted  as 
(a),(b),(c)  and  (d),  are  considered  as  nncorrdated  normal  random  variables.  They  have  a  distri- 
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bution  such  that  the  mean  value  of  the  fixity  factors  is  0.70.  Two  different  c.o.v.,  0.10  and  0.20, 
will  be  used  to  describe  the  dispersion  in  the  values  of  the  stiffiiess  coefficients. 

The  probability  density  function  (PDF)  corresponding  to  the  stiffness  coefficient  of  one  of 
the  flexible  joint  for  three  c.o.v.  is  shown  in  Figure  (1.3).  The  multivariate  density  function 
corresponding  to  the  four  random  stiffness  coefficients  is  obtained  as  the  product  of  the  individual 
PDF  since  the  random  variables  are  assumed  to  be  uncorrelated.  Figure  (1.4)  shows  the  PDF 
for  the  case  of  two  random  stiffness  coefficients. 

If  the  stiffness  coefficients  have  a  normal  distribution,  the  eigenvalues  of  the  system  will 
also  be  normally  distributed.  The  PDF  of  the  first  eigenvalue  is  shown  in  Figure  (1.5).  The 
mean  value  and  the  standard  deviaticm  to  plot  the  curve  were  calculated  with  the  SOPM.  The 
eigenvalues  in  the  horizontal  axis  were  divided  by  the  correspondinsr  eigenvalue  obtained  from 
the  deterministic  eigenproblem  using  the  mean  value  of  the  stifihess  coefficients.  The  c.o.v.  for 
the  random  eigenvadue  are  0.0078  and  0.0158  when  the  c.o.v.  of  the  stiffness  coefficients  are  0.10 
and  0.20,  respectively.  This  indicates  that,  as  in  the  previous  example,  the  structure  is  not  very 
sensitive  to  the  uncertainties  in  the  stiffness  of  the  non-rigid  joints.  Figures  (1.6)  to  (1.8)  show 
the  PDF  for  the  2nd,  3rd  and  4th  eigenvalues.  The  3rd  eigmvaJue  is  the  most  sensitive  to  the 
raundom  variability  of  the  connections’  stiffness.  For  c.o.v.  of  the  stiffiiess  equal  to  0.10  amd  0.20, 
the  respective  c.o.v.  of  the  3rd  eigenvalue  are  0.0146  amd  0.0296.  It  can  be  noted  that  the  2nd 
amd  4th  eigenvalues  do  not  have  a  significamt  dispersion  meaning  that  they  are  not  much  affected 
by  the  uncertainties  in  the  stiffiiess  of  the  connections  at  the  first  floor. 

The  PDFs  for  selected  elements  of  the  lower  four  eigenvectors  are  displayed  in  Figures  (1.9- 
12).  The  results  for  the  1st  eigenvector  atfe  shown  in  Figure  (1.9).  The  ramdom  connections 
have  a  more  pronounced  effect  in  the  modad  dof  associated  with  the  vertical  displacement  of  the 
node  in  the  beaun  c  -  d.  The  two  PDF’s  for  this  dof  depicted  in  the  figure  have  c.o.v.  equal  to 
0.3072  amd  0.6172  for  c.o.v.  of  the  stiffiiess  equal  0.10  and  0.20,  respectively.  Although  the  effect 
is  not  as  pronounced  as  in  the  previous  case,  the  uncertainties  in  the  connections  also  have  an 
important  influence  in  the  rotation  of  joint  (d).  The  PDF’s  associated  with  this  dof  plotted  in 
the  figure  have  c.o.v.  equak^  0.0253  and  0.0517.  On  the  contrary,  the  horizontal  displau:ement 
of  node  (d)  is  relatively  insensitive  to  the  ramdom  connections,  as  evidenced  by  the  sharp  PDF 
for  this  dof  shown  in  the  figure. 

Figure  (1.10)  show  similar  results  but  this  time  for  the  2nd  eigenvector.  The  PDF’s  axe  similar 
to  those  for  the  1st  eigenvector.  The  PDF’s  for  the  modal  displacements  corresponding  to  the  3rd 
eigenvalue  are  presented  in  Figure  (1.11).  In  this  ca^e,  the  largest  dispersion  is  associated  with 
the  horizontal  displacement  of  node  (d).  The  c.o.v.  axe  0.0239  and  0.0485  when  the  respective 
coefficients  for  the  stiffiiess  are  0.10  and  0.20.  Finally,  Figure  (1.12)  displays  the  results  for  the 
4th  eigenvector.  Here  the  vertical  displacement  of  the  interior  node  on  bearnn  (c)-(d)  is  the  modal 
dof  most  aiffected  by  the  ramdomness  of  the  connections.  The  distributions  have  c.o.v.  equal 
to  0.0712  amd  0.1425  when  the  distribution  of  the  stiffness  have  c.o.v.  equal  to  0.10  and  0.20, 
respectively. 
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Table  1.1:  Eigenvalue  Statistics:  Mean  Values. 


Table  1.3:  Eigenvalue  Statistics:  Coefficients  of  Variation. 
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Table  1.4:  First  Eigenvector  Statistics:  Mean  Values. 
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Table  1.5:  First  Eigenvector  Statistics:  Standard  Deviations. 
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Table  1.6:  First  Eigenvector  Statistics:  Coefficients  of  Variation. 
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Fignre  1.4:  Bivariate  Density  Function  of  Stiffness  Coefficients. 
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Figure  1.5:  Probability  Density  Function  of  First  Eigenvalue. 
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Figure  1.6:  Probability  Density  Function  of  Second  Eigenvalue. 
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Figure  1.7:  Probability  Density  Function  of  Third  Eigenvalue. 


Figure  1.8:  Probability  Density  Function  of  Fourth  Eigenvalue. 
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Figure  1.9:  Probability  Density  Fonction  of  Selected  Dof’s  of  First  Eigenvector. 


Figure  1.10:  Probability  Density  f\uiction  of  Selected  DoFs  of  Second  Eigenvector. 


Figure  1.11:  Probability  Density  FVinction  of  Selected  DoTs  of  Third  Eigenvector. 


27 


Chapter  2 

Response  of  Structures  with  Random  Parameters 


2.1  Introduction 

In  the  previons  chapter  we  presented  the  application  of  the  stochastic  finite  element  method  to 
determine  the  statistics  of  the  natural  frequencies  and  modes  of  structures  with  flexible  random 
connections.  In  this  chapter  we  will  apply  the  same  technique  to  calculate  the  statistics  of  the 
response  of  this  type  of  structures  when  they  are  subjected  to  deterministic  dynamic  excitations. 
The  level  of  complexity  of  the  analysis  is  significantly  increased  due  to  the  fact  that  the  perturba¬ 
tion  expansion  that  was  used  in  the  previous  chapter  may  lead  to  a  scdution  that  is  ”  non-uniform!^ 
or  unbounded  in  time.  The  perturbation  expansion  previously  used  is  usually  referred  to  as  the 
” straightforward  expansion^.  The  approodmate  solution  generated  by  this  method  is  only  ap¬ 
plicable  in  the  initial  instants  of  time  due  to  unbounded  terms  appearing  in  the  higher  order 
terms  of  the  expansion.  To  overcome  the  difilculties  associated  with  this  technique,  a  stdution 
methodology  based  on  the  ^method  of  muHtpU  scales”  will  be  devdoped. 

First  we  will  examine  the  deterministic  and  stochastic  response  of  single  dof  systems  in  which 
the  natural  frequency  and  the  excitation  are  functions  of  a  set  of  parameters.  The  parameters 
will  be  conudered  as  detem^stic  perturbations  on  a  set  of  initial  values  in  the  first  case,  and  as 
random  variables  in  the  secmid  case.  This  study  will  not  only  allow  ns  to  comprehend  the  cause 
and  significance  of  the  unbounded  terms  but  the  results  are  important  for  the  final  goal,  i.e,  the 
response  statistics  of  multi-dof  systems  with  random  time-invariant  parameters.  The  latter  study 
is  undertaken  in  the  last  sections  of  this  ch^ter.  For  multi-dof  systems,  the  random  parameters 
will  be  identified  as  the  stiffaess  coefficients  of  the  rotational  springs  representing  the  flexible 
connectimis.  The  case  of  a  sin^  and  multi-dof  system  subjected  to  a  simple  dynamic  loading 
will  be  examined  in  detail  to  illustrate  the  implication  of  the  straightforward  expansion  and  the 
multiple  scales  approach. 


2.2  Single  Degree-of-Freedom  Syetems 

Consider  a  sin^  dof  oscillator  with  natural  frequency  u;,  and  damping  ratio  (  subjected  to  an 
excitation  /(<): 
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^(0  +  2  f  w  iKO  +  »K0  =  fit)  (2.1) 

We  will  assume  that  the  natural  frequency  of  the  oscillator  and  the  excitation  are  function 
of  R  variables  01,03,  ...,aA.  The  reason  for  assuming  this  dependence  will  become  apparent 
later  when  we  calculate  the  statistics  of  the  response.  To  shorten  the  notatimi  we  will  group  the 
parameters  Om  in  a  single  vector  a.  These  parameters  represent  perturbations  in  the  original 
values  of  u  and  /(t).  Without  loss  of  generality,  the  initial  or  unperturbed  state  of  u  and  f{t) 
will  be  associated  with  the  zero  value  of  the  parameters  Om  .  The  response  of  the  oscillator  will 
also  be  a  function  of  the  parameters  Om  mid  equation  (  2.1)  can  be  written  as: 

tj(a, t)  +  2(u;^a,t)  +  rf(a, t)  =  /(o, t)  (2.2) 


2.2.1  The  Straightforward  Expansion 

The  objective  of  the  analysis  presented  in  this  and  in  the  following  sections  is  to  obtain  the  re¬ 
sponse  q(a,  t)  without  having  to  solve  equation  (  2.2)  and  assuming  that  we  know  the  unperturbed 
values  of  a;  and  f(t).  A  straightforward  perturbation  expansion  can  be  used  for  this  purpose.  The 
method  of  straightforward  expansion  is  similar  to  the  second  order  perturbation  method  used  in 
Chapter  4  to  calculate  the  statistics  of  the  eigenvalues  and  dgenvectors.  It  consists  in  obtaining 
an  approximate  expression  for  the  response  by  meaM  of  a  Thylor  series  expansion  in  terms  of  the 
parameters  Am.  Retaining  only  up  to  second  order  terms  in  the  series  we  can  write: 


R  R  ti 

niOyt)  =  »k)(0  +  vLit)  Om  +  I  2  £  »4«(0  Om  On  (2.3) 

msl  msl  nsl 

where  rfo{t)  is  the  response  of  the  system  when  the  parameters  assume  the  initial  or  unperturbed 
values.  The  superscripts  I  and  II  represent,  respectively,  the  first  and  second  order  rates  of 
change  of  the  function  q(a,t)  with  respect  to  the  parameters  Om  evaluated  at  the  initial  state  of 
the  system.  That  is: 


nLi*)  = 


da„ 


QtsO 


vL'M  = 


OsO 


(2.4) 


(2.5) 


dam  dan  I 

The  derivatives  of  the  response  17(0,  t)  with  respect  to  the  parameters  Am  are  not  known  at  this 
stage.  The  goal  of  the  scdution  process  is  precisely  to  obtain  these  derivatives. 


To  simplify  the  notation  in  the  subsequent  developments  the  following  notation  will  be  intro¬ 
duced:  ^ 

mi*)  =  nLi*)<>m  (2.6) 

msl 


H  I3’»«n(0OmOn 


(2.7) 


msl  nsl 
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When  using  perturbation  methods  it  is  customary  to  introduce  a  coefficient  e  to  identify  the 
order  of  the  different  terms  involved  in  the  expansions.  Emplc^ng  equations  (  2.6  )  and  (  2.7  ) 
and  the  perturbation  parameter  £,  the  expansion  (  2.3)  can  then  be  written  as  follows: 

=  Vo(t)  +  e*?i(0  +  (2.8) 

As  it  will  become  evident  later,  for  the  cases  we  are  interested  in,  the  natural  frequency  and 
the  excitation  in  equation  (  2.2)  are  nonlinear  functions  of  the  parameters  Om  •  Ebccept  for  trivial 
cases,  it  is  not  possible  to  obtain  the  exact  closed  form  functional  relationship.  However,  the 
frequency  and  the  excitation  can  also  be  expanded  in  Taylor  series  around  the  initial  values  of 
the  parameters  Om: 


R  R 

U(a)  =  +  E‘*'mnamOn 

mal  msl  nsl 

R  R  R 

/(«,«)  =/o(()  +  E/iWom  +  5  i: 

mal  msl  nsl 

where  the  first  and  second  order  rates  of  change  oi  ut  and  /(f)  are  defined  as: 

^  .  ft(t\  - 

^  f  /mV*/  “  ^ 

0nO  a=Q 


"*  dam 


u;"  * 


d*u; 


da™  da. 


la>0 


.  ,H  _  y/«) 


damdOn 


Introducing  the  following  notation  the  first  and  second  order  correction  terms: 


^  H 

mml 
R  R 


/i(«)  = 

mml 
R  R 


‘^  =  5  E  ;  /l(0  =  5  ^  ^ /mn  “m  « 


mml  nml 


msl  nsl 


(2.9) 

(2.10) 

(2.11) 

(2.12) 

(2.13) 

(2.14) 


and  using  equations  (  2.13-  2.14)  along  with  the  perturbation  parameter  e  the  expansions  in 
equations  (  2.9-  2.10)  can  be  written  in  the  simpler  form: 

u;(o)  =  u>  +  eu;j+e’fa;2  (2.15) 


/(a,f)  =  /o(f)  +  e  hit)  -I-  £»  hit)  (2.16) 

Substituting  equations  (2.8)  and  (2.9-2.10)  in  (2.2)  yields: 

(fkit)  +  emit)  +  e^Ht))  + 
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+  24  (fa>  +  fiwi  +  c'wj)  (jto(l)  +  ernit)  +  e*»?3(0)  + 

+  (w  +  ewi  +  e*4j3)’(i^<)  +  £rfi{t)  +  e*»72(0)  = 

=  /o(0  +  €/i(0  +  eV2(0  (2.17) 

Collecting  terms  of  the  same  order  in  the  perturbation  parameter,  we  obtain  the  following 
hierarchical  equations: 

0(e»)  t 

0(«‘) : 

0(£») : 

where: 

/i(<)  =  Mt)  -  i  xii^noit)  -  2Qunrio(t)  (2.21) 

Mt)  -  flit)  -  24u;3%(0  -  (2ww3  +  «i)  *]to(0  -  24a;i  %(t)  -  2wfa;i  qi(tX2.22) 

It  is  important  to  note  that  the  excitation  terms  fi{i)  for  the  first  and  second  order  equations 
contain,  respectivdy,  the  solutiaiis  of  the  zero  and  first  order  equations.  Since  the  differential 
operator  in  the  three  equations  (  2.18-  2.20)  is  the  same,  these  excitation  terms  associated  with 
the  solutions  of  the  lower  qrder  equations  are  a  potential  source  of  complications.  Indeed,  it  is 
possible  that  resonance  phdh<»neita  may  occur  in  the  higher  order  equations  since  the  forcing 
functions  for  these  equations  are  the  solutions  of  the  equations  of  lower  order.  In  this  case  the 
solutions  will  have  some  terms  that  grow  with  time  which  are  characteristic  of  resonant  responses. 
In  the  tennincdogy  of  perturbation  methods  these  terms  are  called  secular  terms.  As  we  discusse 
later,  this  problem  will  take  place  for  almost  all  excitations. 

To  examine  the  phenomenon  we  proceed  to  study  the  solution  of  the  hierarchical  equations 
(  2.18-  2.20).  Some  of  the  material  in  the  sequel  is  wdl  known  from  the  elementary  theory  of 
ordinary  differential  equations.  However,  it  is  believed  that  its  inclusion  here  is  warranted  because 
it  is  crucial  for  our  discussion  on  the  appearance  of  secular  terms. 


Voit)  +  24<S^t)  +  U^Tj^t)  =  /o(t) 

(2.18) 

mit)  +  24wm(<)  +  w*»/i(t)  =  /i(t) 

(2.19) 

viit)  +  2  4  <*»%(<)  +  <j*i?8(0  =  Mt) 

(2.20) 
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2.2. 1.1  Solution  of  the  Zeroth  Order  Equation 
We  begin  by  examining  the  zero  order  equation  (  2.18): 

Vo(t)  +  =  /o(0 

There  are  two  ways  to  express  the  total  solution  of  equation  (  2.18).  One  way  is  to  write 
down  the  soluticm  as  the  sum  of  an  homogeneous  or  complementary  solutimi  yj^(t)  and  an  inho¬ 
mogeneous  or  particular  solution 


Vo(*)  =  Voit)  + 

The  homogeneous  solution  has  the  form: 

(2.23) 

r^(t)  =  a  Sin(ijdi  +  9) 

(2.24) 

where  a  and  9  are  coefficients  to  be  determined  later  and  Qj  is  the  damped  natural  frequency: 

Qd  =  (2.25) 

The  inhomogeneous  solution  depends  on  the  particular  form  of  the  right  hand  side  term 
in  equation  (  2.18).  The  solution  can  be  obtained  using,  for  example,  the  method  of 

undetermined  coefficients.  The  complete  solutim  can  then  be  written  as: 

ijd(0  =  ae-^^  Sin(Qdt  +  «)  +  (2.26) 

The  coefficients  a  and  9  must  be  obtained  by  applying  the  initial  conditions  to  the  total  response 
given  by  equation  (  2.26). 

Alternatively,  the  solution  of  equation  (  2.18)  can  be  written  in  terms  of  the  response  to  initial 
conditions  and  the  forced  response  r^{t): 

^  (2.27) 

The  response  to  initial  conditions  is  identical  to  equation  (  2.24): 

»j^(t)  =  oo  Sin{Qdt  -I-  flo)  (2.28) 

with  the  difference  that  the  coefficients  oq  and  9o  have  to  be  obtained  by  applying  the  initial 
conditions  only  to  the  solution  defined  by  equation  ( 2.28).  The  forced  response  (t)  is  expressed 
in  terms  of  a  convcdution  integral: 

70^(0=  /  /o(r)h(f-r)dT  (2.29) 

JO 

where  h(t)  is  the  unit  impulse  response  function  which  for  an  underdamped  system  is  defined  as: 

h(t)  =  —e~<^sin((^dt)  (2.30) 
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Therefore,  the  total  solution  becomes: 


%(t)  =  Sin  (Qit  +  tfo)  +  —  /  fair)  Sin  {ij)d{t  -  t))  dr  (2.31) 

Jo 

Based  on  the  analysis  of  the  results  obtadned  so  far,  the  following  observations  about  the  two 
solution  methods  and  the  potential  difficulties  associated  with  them  can  be  made: 

When  the  solution  is  expressed  in  the  form  of  equations  (  2.26)  or  (  2.31),  the  terms  q^(t) 
or  are  harmonic  functions  with  frequency  equal  to  the  natural  frequency  of  the  oscillator. 

Clearly,  these  terms  will  induce  secular  terms  in  the  solution  of  the  first  order  equations. 

When  the  excitation  is  defined  by  a  relatively  simple  expression,  the  response  can  be  obtained 
using  the  method  of  equation  (  2.26).  In  most  situations,  there  will  not  be  any  harmonic  compo¬ 
nents  with  the  oscillator’s  frequency  included  in  the  inhomogeneous  solution  r^^{t).  Therefore, 
in  this  case  the  only  source  of  secular  terms  will  be  the  homogeneous  solution  However, 

one  must  take  into  account  that  both  soluti<ms  become  interconnected  when  the  total  soluticm  is 
used  to  satisfy  the  initial  conditions. 

When  the  integral  in  equation  (  2.31)  is  solved  for  a  specific  loading  function,  the  resulting 
expression  will  have  harmonic  components  with  frequency  equal  to  the  oscillattv’s  frequency.  H 
the  initial  conditions  are  zero,  this  is  the  only  term  that  wiU  introduce  secular  terms. 

2.2.1.2  Solution  of  the  First  Order  Equations 

To  solve  the  first  order  equations  we  substitute  equations  (  2.6)  and  (  2.13)  in  equation  (  2.19), 
and  equate  those  terms  in  both  sides  that  are  multiplied  by  the  same  variable  a^-  The  m-th 
equation  of  first  order  can  then  be  written  as  follows: 

i^it)  +  2^ui^(0  +  =  Hit)  (2.32) 

where:  ^ 

Hit)  -  Hit)  -  2^u;'  ^<)  -  2u>u^Tjoit)  (2.33) 


Substituting  in  the  above  expression  the  zeroth  order  solution  defined  in  equation  (  2.26)  we 
obtain: 

Hit)  =  -2wa;',oe“<®*  ^(l  -  5in(wrft  +  tf)  + 

+  /i(<)  -  2  <4  ({ >iS"(f)  +  "  (2-34) 


The  solution  of  equation  (  2.32)  can  be  written  as  the  sum  of  the  response  due  to  the  initial 
conditions  and  the  forced  response.  For  zero  initial  conditions,  the  latter  term  will  manifest 
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the  limitations  of  the  method.  Therefore,  we  will  focus  our  attention  on  this  term.  The  forced 
response  can  be  obtained  substituting  in  the  Duhamel  integral  the  excitation  defined  in  equation 
(  2.34): 


»7m(0^  =  -  2  w  wjL  a  (l  -  Sin{u;dT  +  $)  h{t  -T)dT- 

-  2  w  a  Co»{ijiT  +  9)  h{t  -  r)  dr  + 

+  /mW  -  +  w»>o"(r))  h{t-T)dT 

Solving  the  first  two  integrals  the  forced  response  becomes: 

\pu)dtCo3{udt  +  9)  -  ru>dt  Sin{Qdt  +  9)  - 


(2.35) 


-  Sin(Qdt)  (p  Co3(9)  -  r  5’in(tf))]  -|- 

+  fLi-r)  -  2  wi,  ri^ir)  +  u  »)^"(r))  h{t  -  t)  dr 

(2.36) 

where  the  coefficients  p  and  r  are: 

p=i-e 

(2.37) 

'~7T^ 

(2.38) 

Note  that  if  damping  is  Neglected,  the  first  two  terms  inside  the  brackets  in 
become  unbounded  due  to  the  terms  linear  in  the  variable  t. 

equation  (  2.36) 

2.2.1.3  Determiniitic  Straightforward  Expansion  for  the  Undamped  Oscillator  Sub¬ 
jected  to  a  Step  Force 

In  order  to  gain  insight  into  the  peculiarities  of  the  approximate  solnticms  obtained  via  perturba¬ 
tion  techniques,  the  specific  case  of  an  undamped  oscillator  subjected  to  a  step  load  of  magnitude 
Fo  will  be  considered.  The  oscillator  is  assumed  to  be  stationary  when  the  load  is  applied.  The 
equation  of  motion  assumes  the  simple  form: 

^at,t)  -H  u;(o)»q(o,l)  =  Fo  (2.39) 

To  express  the  solution  by  means  of  the  perturbation  expansion  in  equation  (  2.8)  we  need 
to  solve  the  hierarchy  of  equations  (  2.18-  2.20).  To  express  the  solution  of  the  zeroth  order 
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equation  (  2.18)  in  the  form  of  equation  (  2.26)  we  need  to  define  the  particular  solution 
By  inspection  of  equation  (  2.39),  it  follows  that: 


The  total  solution  assumes  the  form: 


(2.40) 


=  aSin{Qi  +  6)  +  ^  (2.41) 

and  after  satisfying  the  initial  conditions  it  becomes: 

%(<)  =  §(1-C<w(a;t))  (2.42) 

Obviously,  had  we  used  the  alternative  approach  of  equation  (  2.31)  we  would  have  arrived 
at  the  same  result. 


Next  we  need  to  solve  the  first  order  equation  (  2.19).  We  begin  by  substituting  equations 
(  2.6)  and  (  2.13)  in  equation  (  2.19)  and  equating  the  coefficients  of  the  variables  •  In  this 
way  we  obtain  R  equations,  with  a  typical  equation  having  the  form: 

^(0  =  /m(0  (2.43) 

The  excitaticm  term  is  defined  as: 

/i(<)  =  Swi  f  -  1)  (2-44) 

To  obtain  the  solution  of  first  order  we  will  use  this  time  the  convolutim  integral  approach. 
Since  the  initial  conditions  are  zero,  the  total  first  order  response  is: 

^  /‘(C<»(i)()  -  1)  Mt-T)ir  (2.45) 

^  w  Jo 

Carrying  out  the  integration  leads  to: 

»7m(0  =  ~  {Co${Qt)  +  utSin{Qt)  -  1)  (2.46) 


As  it  was  anticipated,  equation  (  2.46)  exhibits  an  unbounded  behavior  in  time.  Obviously,  this 
is  not  in  conformity  with  the  physical  reality:  the  response  of  an  undamped  oscillator  subjected 
to  a  suddenly  applied  constant  force  should  oscillate  about  the  static  response. 

To  calculate  the  second  order  solution  we  need  to  use  the  zeroth  and  first  order  solution  to 
obtain  each  of  the  terms  qmn(0  second  order  equations.  However,  if  the  solution  ^ven  by 
equation  (  2.46)  is  used,  the  soluticms  of  the  second  order  equations  will  contain  terms  quadratic 
in  t.  In  general,  it  is  possible  to  show  that  if  we  seek  a  solution  based  on  an  expansion  of  order 
n,  the  n-th  term  in  the  expansion  will  have  a  secular  terms  of  the  form  [8]. 
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The  reason  why  the  straightforward  expansion  fails  to  approximate  the  actual  response  of  the 
system  becomes  apparent  if  we  compare  the  total  perturbati<m  solution  with  the  actual  solutitm. 
Using  equations  (  2.42),  (  2.46),  (  2.13)  and  (  2.8)  the  approximate  soluticm  is: 

Fo  ^  jf 

*K0  -  ^  (1  -  Cos(Qt))  +  5^  2t.;4  ^  (Cos(w<)  +  ut  Sin{ui)  -  1)  a„  (2.47) 
whereas  the  exact  soluticm  is  : 

-l(<)  =  (I  -  CmM#)!))  (2.48) 

The  approximate  solution  (  2.47)  tries  to  reproduce  the  exact  solution  with  harmonic  terms 
that  have  the  same  undamped  frequency  Q  of  the  unperturbed  system.  Therefore,  the  approx¬ 
imate  and  exact  solutions  differ  both  in  amplitude  and  phase.  The  phase  difference  makes  the 
two  solution  stand  apart  more  and  more  with  time.  After  a  certain  time,  the  correction  terms  are 
unable  to  correct  the  unperturbed  term.  In  fact,  due  to  the  presence  of  secular  terms  in  the  first 
order  correction  in  equation  (  2.47),  the  fundamental  assumption  of  the  perturbation  method  is 
no  longer  valid. 

2.2.1.4  Response  Statistics  of  the  Undamped  Oscillator  Subjected  to  a  Step  Force 

When  the  response  statistics  are  calculated  with  the  straightforward  expansion,  the  limitations 
of  the  method  will  also  emerge  here,  as  we  will  demonstrate  next.  We  will  consider  again  the 
undamped  oscillator  subjected  to  a  step  loading  function,  but  now  the  parameters  are  zero 
mean  random  variables. 

First  we  will  calculate  the  expected  value  of  the  response  q(a,0-  ^  pointed  out 

before,  the  exact  closed  foon  relationship  between  the  response  T}(o,f)  and  the  parameters  a  is 
not  available.  Hence,  in  ordqp  to  apply  the  expected  value  operator  to  the  response  it  is  convenient 
to  express  y^o,  ()  as  an  expansion  in  terms  of  the  random  variables  Om  •  Retaining  up  to  first 
order  terms  we  have: 


q(«,<)  =  rioit)  +  5^  (2.49) 

msl 

Applying  the  expected  value  operator  £{...}  to  this  expression  we  obtain: 

EiniaM  =  voit)  (2.50) 

Therefore,  the  expected  value  of  the  response  calculated  with  a  first  order  straightforward 
expansion  coincides  with  the  response  of  the  unperturbed  system.  In  other  words,  the  mean 
function  is  equal  to  the  response  of  the  deterministic  system  in  which  the  random  parameters 
Om  are  evaluated  at  thdr  mean  values.  If  instead  of  equation  (  2.47),  we  use  a  second  order 
expansion,  we  would  obtain  second  order  terms  in  equation  (  2.50).  However,  these  terms  would 
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have  quadratic  secular  terms.  For  undamped  systems  the  second  order  terms  grow  indefinitely 
and  soon  the  expansion  becomes  invalid.  In  any  case,  the  result  obtained  with  the  first  order 
expansion,  equation  (  2.50),  is  also  wrong  as  it  is  explained  next. 

To  demonstrate  the  fallacy  in  the  solutimis  obtmned  with  a  straightforward  expansion,  we 
propose  to  carry  out  the  following  conceptual  experiment.  Consider  a  sufficiently  large  number 
of  undamped  oscillator  with  quiescent  initial  conditions.  Let  the  natural  frequencies  of  the 
oscillators  be  random  variables  with  the  same  mean  values.  Suppose  that  the  same  deterministic 
loading  of  limited  duration  is  applied  to  the  oscillators  at  the  same  initial  instant  of  time.  Due 
to  the  different  perturbations  in  the  mean  values  of  the  natural  frequencies  of  the  oscillators, 
their  responses  will  become  more  and  more  out  of  phase  with  time.  Consider  now  an  instant  of 
time  t'  suffidently  far  away  from  the  initial  instant  and  calculate  the  instantaneous  mean  value 
of  the  response  at  this  instant.  The  initially  small  differences  in  the  instantaneous  values  of  the 
individual  responses  grew  so  much  that  at  time  V  they  are  completely  random.  It  is  logical  to 
conclude  that  the  sample  mean  of  the  oscillators'  responses  at  this  time  approaches  zero.  In  the 
limit,  as  we  increase  the  number  of  oscillators  in  the  experiment,  the  sample  mean  estimator  tends 
to  the  mean  or  expected  value  of  the  response  £{q(t')}.  Evidently,  the  mean  value  calculated 
with  a  straightforward  expansion  of  any  order  would  never  be  equal  to  zero  for  the  undamped 
case. 

If  we  repeat  the  experiment  with  a  set  of  damped  oscillators  the  discrepancy  in  the  straight¬ 
forward  expansion  results  described  before  will  not  manifest  with  the  same  intensity.  In  this  case 
the  individual  responses  of  each  of  the  oscillators  would  independently  tend  to  zero  by  virtue  of 
the  damping  in  the  system.  Therefore,  if  we  calculate  the  sample  mean  at  increasing  instants  of 
time  we  will  obtain  that  a  sample  mean  that  decreases  with  time.  If  we  use  the  straightforward 
expansion  to  calculate  the  expected  value  of  the  response,  it  will  yield  a  value  that  does  tend  to 
zero.  However,  it  should  be  clear  that  this  result  is  due  to  the  presence  of  damping  in  the  system 
and  is  not  by  way  of  probabilistic  reasons  as  it  should  be. 

Let  us  turn  back  our  attention  to  the  undamped  oscillator  subjected  to  a  step  load.  We  will 
attempt  to  obtain  the  variance  of  the  response  using  a  first  order  (straightforward)  perturbation 
expansion.  The  variance  of  the  response  can  be  obtained  from: 

=  BMa,  <)  -  E{fHa, !)))»)  (2.51) 

Substitution  of  equation  (  2.49)  in  (  2.51)  yields: 

=  4  S  (Cm(m)  +  ut  Sin<<5()  -  1)’  f;  f;  uW  E{a™<>.}  (2.52) 

msl  n»l 

Taking  into  account  equation  (  2.46)  the  variance  becomes: 

*^5(0  ~  iCo8{Qt)  QtSin{Qt)  -  1)’  53  (2.53) 

^  m»l 
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If  all  the  random  variables  Om  assumed  to  be  uncorrelated,  this  expression  reduces  to: 

^n(0  =  (C<M(wt)  +  a;tS»n(d;t)  -  1)^  <^m  (2-54) 

Examining  the  above  expression  we  conclude  that  the  use  of  a  first  order  straightforward 
expansion  to  calculate  the  response  variance  gives  rise  to  a  quadratic  secular  term. 

2.2. 1.5  The  Method  of  Multiple  Scales 

We  have  shown  that  using  an  approximate  solution  that  does  not  take  into  account  the  change 
in  the  frequency  due  to  the  perturbations  in  the  parameters  of  the  oscillator  gives  rise  to  a 
series  of  problems  that  also  affect  the  response  statistics  calculations.  We  will  examine  another 
approximate  method  that  can  explicitly  take  into  account  the  modificaticm  of  the  frequency.  This 
te^nique,  known  as  the  method  of  multiple  scales,  was  originally  conceived  by  Nayfeh  [8]  and  it 
has  been  successfully  used  in  the  solution  of  weak  nonlinear  systems. 

We  will  consider  again  equation  (  2.2)  for  a  viscously  damped  oscillator  in  which  the  natural 
frequency  is  a  function  of  R  parameters  01,02,... ,ai{.  We  will  again  seek  a  solution  in  the  form 
of  an  expansion  but  now  the  terms  in  the  expansion  will  depend  not  only  in  the  parameters  Om 
but  also  in  products  of  Omlt  etc.  We  can  then  write: 

ij(a,<)  *  y7(oi,O2,...,O|i,oit,O2t,...,OAt,0  (2.55) 

The  solution  can  also  be  written  as: 

q(o,()  —  q(uti«U2,...,Oi{,2b,Ti,T2,  ...,7r)  (2.56) 

where  the  variables  ve  defined  as: 


Tn,  dmt  (2-57) 

« 

We  could  have  considered  that  the  solution  is  also  a  function  of  terms  Omt^*  etc.  This 
is  equivalent  to  keeping  second,  third,  etc.  terms  in  the  expansions.  Since  the  analysis  with  the 
multiple  scales  method  is  more  complicated  than  with  the  strai^tforward  expansion,  in  order  to 
simplify  the  presentation  only  first  order  terms  are  retained  in  the  expansions. 

The  variables  Tq  and  the  set  Ti,  T2,  ...,Ta  can  be  thought  of  as  two  different  time  scales  which 
justify  the  name  of  the  method.  Indeed,  since  Om  ia  a  small  perturbation,  Tm  represents  a  slower 
scale  than  the  zeroth  order  scale  To  =:  t. 

Due  to  the  functional  relationship  assumed  for  the  response  q(Oi0  ia  equation  (  2.56),  the 
time  derivatives  in  equation  (  2.2)  need  to  be  expressed  in  terms  of  the  new  time  scales  To  and 
Ti,T2,  ...Tr  *  Using  the  chain  rule: 


Sill  -  4.  V  ^ 

dl  ~  m  " 


(2.58) 
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(2.59) 


^ 

di^  dTS 


R 

+  2E 
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a*(-) 


®m 


In  order  to  simplify  the  notation  in  the  subsequent  developments  we  introduce  two  differential 
operators: 


y(-) 

dT} 


^  d(...) 

(2.60) 

-  2  f  o 

(2.61) 

Using  the  above  notation  and  the  perturbation  parameter  £  as  a  book-keeping  device  to  keep 
track  the  order  of  the  magnitudes  of  the  different  terms  involved,  equations  (  2.58  -  2.59  )  can  be 
written  as: 


^  ^  a(...) 

dt  ^  dTo  dTi 


(2.62) 


dt»  "  dTS 


(2.63) 


It  is  becoming  clear  now  that  we  will  replace  the  ordinary  differential  equation  of  motion 
(  2.2)  by  a  partial  differential  equation,  thus  complicating  the  original  problem.  However,  this 
apparent  complication  is  far  outweighed  by  the  inherent  advantages  of  the  multiple  scales  method 
[8]  as  we  will  have  occasion  to  verify  it. 


Alike  in  the  straightforward  expansion,  the  natural  frequency  and  the  loading  function  have 
to  be  expanded  in  terms  of  the  variables  a^.  The  series  will  be  truncated  at  the  first  order  and 
we  will  consider  that  the  excitation  depends  only  on  the  lower  time  scale  To  .  Hence,  u;(a)  and 
/(a,  t)  are  approximated  with  the  expansions: 


W(0)  =  W  +  53 


msl 


/(«,<)  =  MTo)  +  2  /iCToK 


msl 


where  the  superscript  I  indicates  the  first  order  rate  of  change  of  the  functions: 


du 

dOm 


;  /to)  =  ^ 


CM) 


(2.64) 

(2.65) 

(2.66) 


Intiodudng  the  fidlowing  definitions: 
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(2.67) 


MTo)  =  E  IUTo)  a„ 


equations  (2.64  -  2.65  )  can  be  written  as: 


(t;(a)  =  u  +  efa;i 


(2.68) 


(2.69) 


/(a, To)  =  /o(ro)  +  £A(ro) 


(2.70) 


The  response  is  also  expressed  as  a  Taylor  series  expansion  truncated  at  the  first  order: 

R 

riia^t)  =  >to(To,Ti,T3,...,Tii)  +  ^  ...,Tr) 

The  rates  of  change  of  the  response  are  defined  in  the  usual  manner: 

vUn.TuTt . r»)  =  ^^^1^ 

Introducing  the  following  notation  for  the  summation  term: 


(2.71) 


(2.72) 


i7i(7btTitT3,...tT|{)  —  •••fTpi)  oin 

mml 


we  can  write: 


»jb(To,ri,r3,..,,rj*)  +  e7i(ro,ri,r3,...,rii) 

Substituting  equations  (  2.62-  2.63),  (  2.69-  2.70)  and  (  2.74)  in  (  2.2)  lead  to: 


+ *«  <® + ‘">)  (^ + '  ^)  (» +*  ■») 


+  (w  +  eu/i)*  (ijb  +  e»7i)  =  /o(To)  +  e  f\{To) 


(2.73) 


(2.74) 


(2.75) 


Collecting  the  terms  of  the  same  order  in  the  perturbation  parameter  e  lead  to  two  hierarchical 
equations  of  motion: 
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where  the  excitation  in  the  first  order  equation  is: 

/i  =  fi(To)  -  -  2ffa;i  -  2a>u;it]b  (2.78) 

Also  in  this  case  one  should  bear  in  mind  that  equation  (  2.77)  actually  implies  Jt  equations  of 
first  order,  each  of  them  associated  with  the  co^cients  in  the  definition  of  rfi  ,  equation 
(  2.73).  Moreover,  it  should  be  noticed  that  here  also  the  excitation  term  in  equation  (  2.78) 
contains  the  response  of  the  equation  of  lower  order,  as  well  as  the  first  two  derivatives  of  the 
response. 

2.2.1.6  Solution  of  the  Zeroth  Order  Equation 

We  will  proceed  to  examine  the  response  of  the  zeroth  order  equation  (  2.76).  As  we  discussed 
earUer,  its  sdution  can  be  written  as  the  sum  of  the  homogeneous  and  inhomogeneous  scdutions 
and  respectively; 

qo(ro,r,,ra,...,rji)  *  q^'(ro,ri,ra,...,rR)  +  ij^^cTo)  (2.79) 

The  homogeneous  scdution  has  the  explidt  form: 

ijg^CTo,  Ti,  Ta,  ...,r#j)  =  a  e'<®^«  5m  (wj  To  +  9)  (2.80) 

It  is  important  to  note  tha^the  coefficients  a  and  0  are  constant  with  respect  to  the  time  scale 
To  .  However,  in  general,  they  do  depend  on  the  scales  Tm  .  That  is: 

a  =  a(Ti,ra,...,TR)  ;  0  =  0(Ti,T2,...,Tii)  (2.81) 

2.2.1.7  Solution  of  tho  First  Order  Equations 

Substituting  in  equatkm  (  2.77)  the  operators  defined  in  equations  (  2.60-  2.61)  and  taking  into 
account  equation  (  2.78)  we  obtain: 


msl 
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Ekiuating  coefficients  of  the  variable  in  both  sides  of  the  equation  leads  to: 

=  (2.83) 

where: 

/m  =  fLC^o)  -  2  ~  ^  -  2w<*;4  (2.84) 

Substituting  the  zeroth  order  solution  defined  by  equations  (  2.79)  and  (  2.80)  in  the  above 
equation  and  carrying  out  the  derivatives  indicated,  equation  (  2.84)  assumes  the  form: 

dTS  ^  dTo  ^  ^ 

=  2u>  +  a  ^  +  (2^*  -  l)u;'  o)  SiniQ^To  +  0) + 

2Q  e’^®^*Cos(u>rfro  +  tf)  + 

+  /i(T,)  -  2{a^.  -  2i>ui  -^"(ro)  (2.85) 

The  first  two  terms  in  the  right  hand  side  of  equation  (  2.85)  will  give  rise  to  secular  terms 
in  the  first  order  sdution.  In  order  to  obtain  a  uniform  expansion,  i.e.  free  of  secular  terms, 
the  method  of  multiple  scales  impose  that  the  coefficients  of  the  trigonometric  functions  with 
frequency  equal  to  the  unperturbed  natural  frequency  must  vanish: 


(2.86) 

d0  da  r 

(2.87) 

Solving  these  equations  we  obtain  the  two  coefficients  a  and  0  that  will  guarantee  that  no  secular 
tenns  will  occur  in  the  expansion,  at  least  up  to  the  first  order. 
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2.2.1.8  Deterministic  Multiple  Scales  Expansion  for  the  Undamped  Oscillator  Sub¬ 
jected  to  a  Step  Force 

Before  proceeding  with  the  calculations  of  the  statistics  of  the  response  it  will  prove  advantageous 
to  apply  the  method  of  multiple  scales  to  calculate  the  deterministic  response  of  an  oscillator  sub¬ 
jected  to  an  excitation  described  by  a  mathematically  simple  expression:  a  step  forcing  function. 
For  simplicity,  damping  will  be  neglected  and  the  oscillator  is  assumed  to  be  initially  at  rest.  The 
equation  of  motion  reduces  to  equation  (  2.39) 

=  Fo 


The  solution  of  order  zero  is: 

i]u(ro,ri,r2,...,r„)  =  aSin{Qt  +  9)  +  ^  (2.88) 

U) 

For  this  particular  case  the  equation  of  first  order  becomes: 

»«  =  2w  (a  ^  Sin{u)To  -I-  8)  -b 

CosiQTo  +  8)  -  2u;'  ^  (2.89) 

aim 

The  problematic  terms  are  the  two  trigonometric  functions  with  frequency  Q.  Therefore,  the 
conditions  for  elimination  of  secular  terms  require  that: 

(MO) 

^  =  0  (2.91) 

FVom  the  second  equation  it  follows  that: 

0  =  00  +  0{e^)  (2.92) 

Note  that  the  coefficient  o  is  constant  with  respect  to  the  time  scales  Tm  •  However,  it  may 
depend  on  time  scales  of  the  second  order  (which  we  did  not  include  in  the  analysis).  This  is 
indicated  in  equation  ( 2.92)  with  the  term  O(e^).  The  other  coefficient  9  can  be  found  considering 
equations  (  2.90)  and  (  2.92).  It  follows  that: 

*=  +  (2.93) 

msl 

The  soluti<ms  of  the  R  equations  of  first  order  consist  of  the  particular  soluticms  only.  Therefore, 
by  inspection  of  equation  (  2.89)  we  can  write: 
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»?m  =  -2«jL§  (2.94) 

Substituting  the  coefficients  a  and  6  from  equations  (  2.92-  2.93)  in  equation  (  2.80)  and 
adding  the  zeroth  and  first  order  solution  according  to  equation  (  2.74),  we  obtain: 

iKO  =  ao5*n[wt  +  S ‘^m  «m  <  +  ^  |  (2.95) 

\  m*l  /  \  m=l  / 

The  final  task  is  to  obtain  the  undetermined  coefficients  oq  and  Evaluating  the  above  expres¬ 
sion  in  17(0)  =  i)(0)  =  0  it  is  easy  to  show  that  the  total  solution  is: 

“  I  -  CoS^Qt  +  (2.96) 

It  is  interesting  to  compare  the  multiple  scales  solution,  equation  (  2.96),  with  the  straightfor¬ 
ward  expansion  solution,  equation  (  2.47).  We  immediately  see  that,  unlike  the  straightforward 
expansion,  the  multiple  scales  solution  does  include  a  first  order  correction  in  the  unperturbed 
natural  frequency  u>. 

2.2.1.9  R4S8poii«e  Statistics  of  the  Undamped  Oscillator  Subjected  to  a  Step  Force 

The  case  in  which  the  parameters  Om  are  random  variables  with  zero  mean  will  be  examined 
next.  The  objective  is  to  calculate  the  mean  value  function  and  the  variance  function  of  the 
response  q(t). 

To  obtain  the  mean  value  we  apply  the  expectation  operator  to  equation  (  2.96): 

cw*)}  =  ®  {§  (‘  - 1  r  "i  “")  (i  -  (*• + £  “"'))} 

This  expression  can  be  re-accommodated  as  frdlows: 

^w')}  =  §  (i  -  ®|(i  - 1  (2-98) 

Before  continuing  it  is  relevant  to  underscore  the  fact  that  even  though  we  used  a  first  order 
expansion,  the  expression  for  the  expected  value  will  differ  from  the  response  of  the  unperturbed 
system.  This  is  in  contrast  with  the  results  obtained  with  the  straightforward  expansion. 

We  need  to  solve  the  multiple  integrals  associated  with  the  following  expected  values: 


=  J  ^  •  J  f  (Ol  02  On)  Cos  Om  dtti  rfttj  ...  dOR  (2.99) 

CoS^t+  = 

=  y^  •••  y^  /  (oi  aj  ...  an)  ^  o«^  cos  *  j  dai  da^ ...  doR 

(2.100) 

where  the  fanGtion  /  (01,03,  ..,an)  is  the  jdnt  probability  density  function  of  the  random  vari¬ 
ables  On  .  The  evaluation  of  these  integrals  is  rather  laborious.  The  linal  results  are: 

£|c<M^u;t-|-  =  H  Cosiut)  (2.101) 

Cos^w*  +  = 

=  -*5»n(u«)  53  ((<^mWm)  e“  *  n  e"  2  j  (2.102) 

msil  \  na:l;n|tm  / 


Therefore,  the  expected  value  of  the  response  can  be  written  as: 


E{fKt))  =  S  (1  -  iCosiCH)  -b  Q{t)Sin{Qt))  R{t)) 


where  the  following  notation  is  used: 


mat 

(}(()  =  I  P< 

w 

R{t)  = 


(2.103) 


(2.104) 


(2.105) 


(2.106) 


Next  we  will  determine  the  variance  of  the  response.  We  begin  by  calculating  first  the  mean 
square  value  of  the  response.  Using  equation  (  2.96),  we  can  write: 
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+  2  -  1  j  Cos  +  53  Wm  am  - 

“  (l  S  ‘^m  “m  -  ij  Cos^  ^U>t  +  53  ‘^m  Om 


(2.107) 


The  double  summatioii  terms  in  the  variables  amOn  have  been  neglected  in  the  above  equation. 

Applying  the  expected  value  operator  on  equation  (2.107)  and  recalling  that  the  variables 
have  zero  mean,  we  obtain: 


emn  = 


(i)’  +  ^^{(1 1“"“”  -  *) 

”  -1^C’0S*^W«+ 


(2.108) 


To  proceed  further  we  need  to  calculate  the  expected  values  of  expressions  that  involve  the  square 
of  the  cosine  function.  This,  in  turn,  requires  to  evaluate  the  following  multiple  integrals: 


E 


|cos* 


A 

wt  +  51 

nwl 


"mOm* 


109) 


=  j  /  (ai  03...  or)  Cos^  fw*  +  53  "rnttm  A  doida2...doj|  (2. 

\  9f|sl  / 


=  ££•••£  /(Oi  03  ...  or)  CoS*  ffciX  +  51  ‘*'m  dOi  da2  ...  dOR 

00  oo  00  Vmsl  /  \  msl  / 


(2.110) 


After  performing  the  evaluation  of  the  integrals  we  obtain: 
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^  n  C03(2w<)^  (2.111) 

^  I  (  IZ  ""*)  +  51  ‘*>m  Om  I  = 

=  -tSin(^t)  ^  I  (<ymW«)*e“^(‘''"""*)  JJ  e-2(<^""i0  |  (2.112) 

msl  \  n=l;n^m  / 

Using  the  above  results  along  with  the  expected  values  of  the  expressions  in  equation  (  2.101- 
2.102)  and  (  2.104-  2.106),  the  mean  square  response  can  be  written  as: 

=  (§)'  (I  -  2(C<»(i>0  +  QWSMat))  .-i'”*’  + 

+  ^(C<w(2w()  +  ^(*)5in(2u)())e-^^‘*)  (2.113) 

We  note  that  as  time  passes,  the  mean  square  response  tends  to  the  value  of  the  unperturbed 
system.  That  is, 

to£{^<)’}  =  |(|)’  (2.114) 

Finally,  to  obtain  the  response  variance  function  we  have  to  substitute  equations  (  2.103)  and 
(  2.113)  in  the  following  expression: 

••  <7^,,  =  £(.^1)’)  -  (£{>K<)))’  (2.115) 

2.3  Multi-Degree-of-EVeedom  Systems 

We  will  consider  structural  systems  in  which  the  mass  and  stiffness  matrices  as  well  as  the 
excitation  vector  are  function  of  the  set  of  parameters  01,02,  ...,aA.  These  variables  represent 
perturbations  on  the  values  of  a  certain  set  of  parameters  of  the  system.  We  are  specially 
interested  in  the  case  in  which  the  parameters  0^3  represent  random  variations  about  the  mean 
values.  Using  a  vector  o  to  collect  the  set  of  random  variables,  the  equations  of  motion  can  be 
written  as: 

[M(«)]  ^a,t)  +  [C]  «(«,<)  -I-  [Kia)]  •(«,<)  =  /•(«,<)  (2.116) 

The  excitation  vector  is  assumed  to  be  a  function  of  the  parameters  o^  in  order  to  include 
excitations  of  seismic  origin  in  the  analysis.  In  this  case  the  vector  F{a,  t)  assumes  the  form: 
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(2.117) 


F{a,t)  =  -  (Af(a)j  ria(0 

where  Xg{t)  is  the  acceleration  of  the  base  and  r  is  the  vector  of  influence  coefficients. 

The  solution  of  the  equations  of  motion  (  2.116)  can  be  obtained  applying  perturbation 
techniques  to  the  coupled  equations  of  motion  or  to  the  uncoupled  system  of  equations.  We  will 
examine  both  cases. 

2.3.1  Perturbation  Method  Applied  to  the  Coupled  Equations  of  Motion 

To  obtain  a  solution  in  the  form  of  a  perturbati<m  expansion  it  is  necessary  to  replace  all  the 
system  matrices  as  well  as  the  response  and  excitation  vectors  by  their  respective  Taylor  series 
expansion  with  respect  to  the  variables  ot.  Keeping  only  flrst  order  terms,  we  can  write: 

[Ar(o)]  =  [K]  +  (2.118) 

msl 

=  [fll  +  f;  [<]  c,„  (2.119) 

msl 

=  !•(<)+  (2-120) 

msl 

1l(O,0  =  1lo(<)  +  (2.121) 

msl 

To  shorten  the  notation  the  following  terms  are  defined: 

(jr.)  =  •£  [Ari]  (2.122) 

^  msl 

[Mil  =  x;  K1  (2‘“) 

msl 

/•iW  =  EfI.  (<)“..  (SIS'!) 

msl 

•iW  =  E  •£.(«)<>»  (2.125) 

msl 

Equations  (  2.118-  2.121)  can  then  be  expressed  as: 

[A(a)J  =  [K]+e  [Ki]  (2.126) 

[M(o)]  =  [M]  -^e  [Ml]  (2.127) 
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Fia,t)  =  m  +  eFiit) 


(2.128) 


ii(a,0  =  «o(0  +  (2.129) 

Substituting  the  expansions  in  equations  (  2.126-  2.129)  in  the  equations  of  motion  (  2.116)  gives: 

(IM)  +  £  [M, ])(«(()  +  »,(())  +  lCl(«o(f)  +  E*i(<))  + 

+  (1^1  +  £  (Jfil)  (lio(l)  +  £«l(l))  =  Hi)  +  £fi(l)  (2-130) 


The  following  hierarchical  equations  are  obtained  after  collecting  terms  of  the  same  order  in  the 
perturbation  parameter  e: 

Oie^)  =  [M]io(0  +  [C]ito{t)  +  [^]«o(0  =  ^(0  (2.131) 

0(£i)  =  (M]ii(0  +  +  [^]«i(0  =  Flit)  (2.132) 

wll6r6! 

Flit)  =  Flit)  -  (M,]io(0  -  [Ki]itoit)  (2.133) 

It  is  important  to  keep  in  mind  that  equation  (  2.133)  actually  implies  R  systems  of  equations. 
This  is  so  because  the  vectors  «(t)  and  Fit)  are  defined  as  linear  combinations  of  R  vectors  u^(t) 
and  Ft^it),  respectively.  Note  also  that  the  two  sets  of  equations  (  2.131-  2.132)  possess  the  same 
matrix  differential  operator  and  only  differ  in  the  excitation  term. 

It  will  be  assumed  that  the  damping  matrix  corresponds  to  a  classical  damping  model  or  that 
the  system  has  Rayldgh’s  damping.  Hence,  the  zeroth  and  higher  order  systems  can  be  decoupled 
using  the  following  transformations: 


Mt)  =  (#]*»){0  (2.134) 

«i(0  =  [♦l*h(0  (2.135) 

The  eigenvectors  in  [4]  are  those  of  the  undamped  unperturbed  system.  They  are  obtained  from 
the  soluticm  of  the  eigenproblem: 


{[K]  -  A.  [M])  =  0 


(2.136) 


Substituting  the  expressions  ( 2.134-  2.135)  in  equations  (  2.131-  2.132)  and  premultiplying  on  the 
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right  by  the  transpose  of  the  matrix  [^]  leads  to  the  following  two  sets  of  uncoupled  differential 
equations: 


Oi£°):  Tjio{t)  +  2(i<Dii^(t)  +  =  iV.o(t)  (2.137) 

:  %(0  +  %(t)  +  uf  i^,(0  =  (2.138) 

where  the  modal  excitation  terms  are  given  by: 

NM  =  (2.139) 

Ni,(t)  =  '  If  flb(0  -  »*^nb(0  (2.140) 

The  vectors  Ft  and  in  the  last  equation  are  defined  as: 

rj  =  (Mil  [♦]  (2.141) 

*f  =  if  |Jf.l  1*1  (2.142) 


Recalling  that  the  vector  «i(f)  is  defined  as  an  expansion  in  terms  of  the  parameters  am  &nd  that 
i^(t)  and  «i(t)  are  related  through  the  transformation  in  equation  (  2.135),  we  can  express  %(t) 
as: 

R 

niit)  =  (2.143) 

msl 

Each  of  the  first  order  terms  identified  with  the  subscript  1  in  equation  (  2.138)  is  defined  as 
a  summation  of  derivatives  Multiplied  by  the  variables  •  Substituting  these  expressions  and 
equating  the  coefficients  of  the  same  variable  Om  one  obtains  R  equations  with  the  form: 

i{.(()  +  2{,i),  iJ„(()  = 

=  E  *n^-(<)  -  r.-  «»(<)  -  E  Tifc.  «^,(<) 

msl 

N 

-  *iim  Vioit)  -  53  *^*o(*) 

in  which  and  r,-A„  are  the  elements  of  the  following  two  vectors; 

If.  =  K1  w 


(2.144) 

(2.145) 
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(2.148) 


*L  =  [<]  [*i 


It  is  seen  that  the  excitation  in  the  i  —  th  equation  of  first  order  includes  the  i  —  th  modal 
acceleration  qio(0  and  displacement  qio(t)  obtained  from  the  solution  of  the  zeroth  order  equation. 


These  functions  will  act  as  resonant  loading  for  the  first  order  equations,  therefore  introducing 
secular  terms  in  the  solutions.  Liu,  Besterfield  and  Belytschko  [39,40]  proposed  to  use  a  procedure 
based  on  a  numerical  Fourier  transform  to  filter  out  the  harmonic  components  in  the  excitation 
of  the  first  order  equation.  Eliminating  in  this  way  the  terms  that  lead  to  secular  terms  they 
were  able  to  obtain  a  uniform  expansion.  It  should  be  pointed  out  that  the  secular  terms  appear 
because  the  solution  is  only  valid  for  short  periods  of  time.  Although  by  eliminating  these  terms, 
the  expansion  is  rendered  uniform,  this  process  does  not  intrinsically  improve  the  approximate 
solution.  Indeed,  the  method  still  tries  to  estimate  the  exact  response  using  terms  with  frequencies 
equal  to  the  natural  fluencies  of  the  unperturbed  system. 

2.3.2  Perturbation  Method  Applied  to  the  Decoupled  Equations  of  Motion 

In  this  section  we  will  present  an  alternative  way  to  apply  the  perturbatimi  technique  for  the 
solution  of  the  equations  of  motion  (  2.116).  Although  we  may  not  know  this  matrix,  it  is 
always  possible  to  decouple  the  equations  of  motion  using  the  eigenvector  matrix  [^(a)]  in  the 
trmisformation: 


a(a,0  -  [♦(«)]  (2.147) 

The  only  conditions  are  those  usual  in  modal  analysis,  namely  that  the  system  must  be  linear  and 
classically  or  proporticmally  damped.  If  the  system  has  Rayleigh  damping,  that  is  if  the  matrix 
[C]  is  defined  as  a  linear  combination  of  [iir(a)]  and  [^(a)] ,  the  equations  of  motim  can  also  be 
decoupled,  but  the  modal  dfimping  ratios  would  be  function  of  the  parameters  am  •  We  are  not 
going  to  study  this  case  and  we  will  concentrate  in  the  case  in  which  the  damping  ratios  are 
all  independent  of  Om  •  The  uncoupled  equations  of  motion  obtained  from  the  transformation 
(  2.147)  are: 


^•(a,0  +  2fju;<(o)^(o,t)  +  Ui{afrii{a,t)  =  Ni{a,t) 
where  the  i-th  generalized  force  is: 


(2.148) 


iV<(a,t)  =  #.(a)^  J(a,0 


(2.149) 


This  method  is  more  amenable  for  analytical  study  using  the  results  of  Section  2.2.2.  Therefore, 
it  is  adopted  to  examine  in  detiul  a  specific  example. 


52 


2.S.2.1  Deterministic  Response  of  Multi  Dof  System  Subjected  to  a  Step  Force 

Consider  an  undamped  multi-dof  system  in  which  the  excitation  is  defined  by  a  mathemati¬ 
cally  simple  expression.  The  same  step  loading  function  previously  used  for  single  dof  system  is 
selected.  The  excitation  vector  is  then: 


Fit)  =  FoSUit)  (2.150) 

where  5  is  a  vector  with  constant  coefficients,  Fq  is  the  magnitude  of  the  load  and  U{t)  is  the 
Heaviside  function.  Ignoring  the  system  damping  the  equation  (  2.148)  reduces  to: 

q.(a,t)  -I-  uiia)^ riiia,t)  =  Niia,t)  (2.151) 

where: 

Niia,t)  =  Foikiia)'^  suit)  (2.152) 


A  crucial  simplification  will  be  made  at  this  point.  We  will  neglect  the  effect  of  the  parameters 
on  the  eigenvectors  ^a)  in  the  definition  of  the  modal  excitation  Ni  .  A  similar  assumption  was 
used  by  Prasthofer  and  Beadle  [11]  in  their  study  of  the  response  of  a  system  with  a  single  random 
parameter.  This  assumption  will  greatly  simplify  the  calculations  of  the  response  statistics  in  the 
next  section.  We  can  then  write: 

Nii<a,t)  a  Niit)  =  Fo^JSUit)  (2.153) 

Recalling  the  expression  (  2.96)  in  Section  2.2.2  for  the  response  of  a  single  dof  system  sub¬ 
jected  to  a  step  force  we  have  that: 

^  ^  S ®"*)  ^  «"*  (2154) 

where  « 

Ni  =Fo^fS 

To  obtain  the  response  in  terms  of  the  physical  coordinates  ii(t)  we  will  make  use  of  the  same 
simplification  as  in  equation  (  2.153): 


N 

«i(0  =  (2-155) 

2.3.2.2  Response  Statutics  of  a  MultUDof  System  Subjected  to  a  Step  Force 

We  will  consider  now  the  problem  of  determining  the  statistics  of  the  response  of  a  multi-dof 
system  in  which  the  mass  and  stiffness  coefficients  are  function  of  random  variables  01,03,  ...,o/{. 
As  before,  these  variables  represent  random  perturbation  about  the  mean  values.  The  mean 
values  will  be  assumed  equal  to  zero. 
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The  mean  value  of  the  displacement  of  the  i~th  dof  can  be  calculated  with  equation  (  2.155): 

/V 

(2.156) 

i=i 

Using  the  results  of  Section  2.2.2  the  expected  value  of  the  modal  displacement  rjj{t)  is  defined 
as: 

=  §  (l  -  (Cos{Qjt)  +  Qj(t)Sin{(jjt))e-2^i^)  (2.157) 

yjj  '  ' 

where: 

R  j 

Py  =£(<'"“!.)  (2158) 

msl 

Qj(t)  =  —  Pit  (2.159) 

To  obtain  the  variance  of  the  response  vector  «(<)  we  proceed  as  follows.  The  variance  of  the 
i  —  th  displacement  can  be  calculated  according  to: 

-  («{«.(<)})’  (2.160) 

Recalling  equation  (  2.156)  the  second  term  in  the  right  hand  side  of  the  previous  equation  can 
be  written  as: 

(£{«((<)»’  =  E  E  *iy  ♦.*  £{>!,(<)}  £{-»(<)}  (2.161) 

jail  isl 

Moreover,  according  to  equ^ion  (  2.155): 

N 

«.(<)’  =  im^ii*.*»7i(07*(0  (2.162) 

kml 

It  is  convenient  to  write  this  expression  as: 

+  2  ^  53  (2.163) 

te>+l 

The  first  term  in  equation  (  2.160)  can  be  obtained  calculating  the  expected  value  of  the  above 
expression: 

^  51  ^ij^ikE{rii{t)nk{t)}  (2.164) 

i»i  j»i  *»i+i 
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The  first  term  of  this  equation  is  the  mean  square  value  of  the  j  -  th  modal  displacement.  This 
can  be  calculated  using  equation  (  2.113)  for  a  single  dof  system: 


=  (I)’  (I- 

-2(Cos(wjt)  +  2Qi{t)Sin{uit))e-hPi^  ^ 

+  i  (Cos(2u>jt)  +  4Qj{t)Sin{2ujt))e~^^^^^  (2.165) 

where  Pj  and  Qj(t)  are  defined  in  equations  (2.158  and  2.159). 

The  terms  in  the  second  summation  in  the  mean  square  value  of  the  response  are  the  expected 
values  of  the  product  After  some  algebraic  manipulations,  the  expression  for  these 

expected  values  reduces  to: 

cWOlW}  =  (|§)  (1- 


+  Bjkit)  e~  i  <*  +  Cjk{t)  e"  i  **  (2. 166) 

where  the  following  definitions  have  been  introduced: 

Ar,(t)  =  C<M(i>rt)  +  *<  (^  +  ^)  ««(<*<)  (2.167) 

Br.W  =  i  (cm((0,  +<5.)()  +  2'  (^  +  ^)  S«((*  +'2.)<))  (2.168) 

C„(#)  =  1  (c<w((ii,+6),)l)  +  21  (^  +  ^)s<»((iv +".)<))  (2.169) 

<J-.  =  (2170) 

Mr,  =  (2.171) 

mm\ 
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(2.172) 


mml 

*-=  E-^  ("J- +"!.)’  (2173) 

mst 

O..  =  f:»i  K-ui.)’  (2.174) 

msl 


Finally,  the  variance  of  the  response  is  calculated  from  equation  (  2.160)  with  the  two  terms 
defined  by  equations  (  2.161)  and  (  2.164)  along  with  equations  (  2.165)  and  (  2.166). 


2.3.2.3  Numeriesd  Examples 

The  plane  frame  shown  in  Figure  (1.1)  will  be  used  to  present  the  calculation  of  time-histwy 
response  statistics  for  a  multi-dof  system. 

We  will  examine  the  undamped  response  to  a  step  loading  function  acting  on  joint  (a)  in 
the  horizontal  direction.  The  stiffiiess  coefficients  of  the  connections  at  the  beam’s  ends  are 
considered  as  nncorrdated  normal  random  variables.  The  mean  value  of  the  stiffiiess  coefficients 
is  such  that  the  corresponding  fixity  factor  is  0.50,  and  two  different  values  of  c.o.v.,  0.10  and 
0.20,  will  be  used  to  describe  the  degree  of  uncertainty  in  the  distribution. 

Figure  (2.1)  shows  the  time  history  for  the  expected  value  of  the  horizontal  displacement  of 
joint  displacement  of  joint  (a).  The  curves  were  normalized  by  dividing  them  by  the  maximum 
value  of  displacement  of  the  system  with  (deterministic)  rigid  joints.  The  time  ordinates  were 
also  divided  by  the  fundamental  natural  period  of  the  rigid  system.  The  effect  of  the  uncertainty 
associated  to  the  connection  stiffisesses  can  be  seen  in  this  figure:  instead  of  oscillating  with 
constant  amplitude,  the  expected  value  of  the  displacement  presents  a  decrement  for  successive 
cycles.  As  time  increases,  the  expected  displacement  tends  to  a  constant  value.  This  effect  can 
be  also  seen  in  Figure  (2.2)  which  shows  the  expected  value  and  the  mean  square  value  of  the 
horizontal  displacement  of  joint  (a),  for  a  c.o.v.  equal  to  0.20.  As  time  increases,  the  mean  square 
value  also  tends  to  a  constant  value. 
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Figure  2.1:  Expected  Value:  Horizontal  Displacement  of  Joint  (a). 


Figure  2.2:  Expected  Value  and  Mean  Square  Value:  Horizontal  Displacement  of  Joint  (a). 
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Chapter  3 

Conclusions  and  Recommendations  for  Future  Work 


3.1  Conclusions 

A  finite  element  model  of  a  beam  element  was  formulated  to  incorporate  the  effects  of  the  flexi¬ 
bility  and  eccentricity  of  the  end  connections  on  the  dynamic  behavior.  A  variational  formulation 
was  employed  to  obtain  the  mass  amd  stiffness  matrices  that  consider  both  effects.  The  system 
matrices  were  expressed  as  the  sum  of  standard  finite  element  matrices  plus  correction  matrices. 


The  model  developed  was  used  in  several  numerical  examples  to  examine  how  the  dynamic 
properties  and  response  of  structural  elements  are  affected  by  the  finite  size  and  flexibility  of 
the  connections.  It  was  observed  that  the  lower  natural  frequencies  diminish  as  the  connection 
stiffness  reduces,  whereas  their  values  increase  with  the  joint  size.  It  was  shown  that  a  structure 
with  initially  rigid  joints  and  non-coincident  natural  frequencies,  may  present  frequency  crossing 
phenomena  when  the  stiffness  of  the  joints  is  changed.  For  structures  subjected  to  localized  har¬ 
monic  excitations  it  was  shown  that  changes  in  connection  stiffness  can  alter  the  frequency  which 
presents  the  maximum  amplification.  The  case  of  structures  subjected  to  seismic  excitations 
was  also  considered.  Severa^design  response  quantities  of  a  ten-story  framework  were  calculated 
using  a  response  spectrum  method  and  a  modified  SRSS  rule  to  combine  the  modal  responses. 
It  was  shown  that  the  joint  characteristics  (stiffoess  and  size)  have  a  more  pronounced  effect  on 
the  design  displacements  than  on  the  design  forces.  For  the  structure  considered,  it  was  observed 
that  an  increase  in  the  size  cf  the  connections  provokes  an  increase  in  the  design  base  shear  and 
a  decrease  in  the  design  displacements. 

Expressions  to  calculate  the  derivatives  of  eigenvalues  and  eigenvectors  with  respect  to  the 
fixity  factors  were  derived  explicitly,  and  they  were  used  to  study  the  eigenvalue  sensitivity  to 
changes  in  the  stiffness  of  the  joints.  It  was  shown  that,  for  a  given  structure  and  depending 
on  which  joints  are  perturbed,  there  is  a  maximum  sensitivity  interval  for  each  eigenvalue.  It 
was  also  observed  that  it  is  possible  to  determine  a  certain  pattern  of  perturbations  that  will 
produce  the  maximum  variation  of  a  given  eigenvalue.  In  general,  this  is  achieved  by  introducing 
on  selected  joints,  perturbations  which  do  not  have  necessarily  the  same  magnitude. 

To  take  into  account  the  uncertainty  in  the  values  of  the  connection  stiffness,  they  were 
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modelled  as  random  variables.  A  second  order  perturbation  approach  was  used  to  formulate  a 
stochastic  finite  element  model,  which  was  used  to  predict  the  statistics  of  the  eigenvalues  and 
eigenvectors.  The  results  were  validated  by  comparing  them  against  a  Monte-Carlo  Simulation. 
An  excellent  agreement  was  observed.  The  numerical  examples  showed  that  the  level  of  dispersion 
induced  by  the  uncertainty  of  the  connection  stiffness  varies  with  the  eigenvalue  number  and  is 
different  for  the  different  components  of  the  eigenvectors.  For  the  example  considered,  it  was 
found  that  an  almost  linear  statistical  relation  exists  between  the  eigenproperties  and  the  stiffness 
of  the  joints. 

The  calculation  of  the  statistics  of  the  time-history  response  of  single  and  multi-dof  structures 
with  random  parameters  was  studied.  It  was  shown  that  the  use  of  the  straightforward  expansion 
leads  to  results  that  are  not  in  agreement  with  the  physical  reality,  from  a  probabilistic  point  of 
view.  A  solution  methodology  based  on  the  method  of  multiple  scales  was  proposed  to  overcome 
the  shortcomings  associated  with  the  straightfcx’ward  expansion  technique.  The  formulation 
developed  was  used  to  determine  the  response  statistics  of  a  structure  subjected  to  a  simple 
loading  case. 

3.2  Recommendations 

It  is  hoped  that  the  work  in  these  two  reports  can  open  the  door  for  future  studies  that  will 
further  enhance  the  understanding  of  the  phenomenon.  Some  of  the  possible  extensions  and 
generalizations  of  this  study  are  the  following: 

The  combined  effect  of  the  connections  flexibility  associated  with  rotational,  torsicmal  and 
extensional  degrees  of  freedom  can  be  included  in  the  finite  element  model. 

The  energy  dissipatimi  in  the  joints  can  be  taken  into  account  by  means  of  rotational  dampers 
at  the  joints.  An  equivalent  viscous  damping  model  can  be  used  to  define  the  coefficients  of 
the  dampers.  The  corresponding  damping  matrix  may  be  developed  following  the  methodology 
presented  in  Chapter  2  of  tlTe  first  report. 

The  sensitivity  analysis  presented  in  Chapter  3  of  the  first  report  can  be  also  be  extended  to 
include  the  effect  of  localized  damping  in  the  joints. 

The  effect  of  the  correlati<m  among  the  random  variables  representing  the  stiffiiess  of  the 
connections  can  be  included  in  the  stochastic  eigenvalue  problem  in  Chapter  1  of  the  second 
report. 

The  methodology  develt^ed  in  Chapter  2  of  the  second  report  can  be  generalized  to  include 
more  complicated  loading  functions  than  those  used  in  this  study.  For  instance,  the  case  of  a 
linearly  varying  excitation  can  be  very  useful  to  study  the  response  due  to  loadings  defined  by  time 
histories  formed  by  linear  segments,  such  as  earthquake  accelerograms.  In  this  case  one  should 
also  consider  the  case  in  which  the  initial  conditions  are  different  from  zero  at  the  beginning  of 
the  time  interval  in  which  the  load  is  divided. 

The  effect  of  damping,  neglected  in  the  analysis  in  Chapter  2  of  Part  II,  can  be  taken  into 
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account  in  the  calculation  of  the  response  statistics.  This  is  specially  important  to  be  able  to  treat 
seismic  excitations  or  other  non-impulsive  loads.  The  modal  damping  ratios  can  be  regarded  as 
deterministic  variables  or,  in  a  more  general  but  also  more  complicated  case,  as  random  variables. 


Finally,  the  nonlinear  characteristics  of  the  moment-rotation  curves  should  be  taken  into 
account.  In  this  case  the  finite  element  matrices  developed  in  Chapter  2  of  Part  I  can  be  used  as 
“tangent  matrices”  for  the  nonlinear  analysis. 
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